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, We consider the proximity effect in multiterminal ferromagnet / superconductor (FSF) hybrid 

structures in which two or three electrodes are connected to a superconductor. We show that 
£SJ \ two competing effects take place in these systems: (i) pair breaking effects due to the response 

to the exchange field induced in the superconductor; (ii) a reduction of the superconducting order 
parameter at the interface that takes place already in NS junctions. We focus on this second effect 
that dominates if the thickness of the S layer is small enough. We consider several single-channel 
electrodes connected to the same site. We calculate the superconducting order parameter and the 
local density of state (LDOS). With two ferromagnetic electrodes connected to a superconductor 
we find that the superconducting order parameter in the ferromagnetic alignment is larger than the 
superconducting order parameter in the antiferromagnetic alignment (Af > Aap), in agreement 
with [Eur. Phys. J. B 25, 373 (2002)]. If a third spin polarized electrode is connected to a 
superconductor we find that Ap — Aap can change sign as the transparency of the third electrode 
■ increases. This can be understood from the fact that the superconducting order parameter is reduced 

if pair correlations among the ferromagnetic electrodes increase. If the two ferromagnetic electrodes 
are within a finite distance we find Friedel oscillations in the Gorkov function but we still obtain 
A F > A AF . 
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I. INTRODUCTION 



The manipulation of entangled states of electrons in condensed matter devices has focussed an important interest 
recently. The ground state of a superconductor is a condensate of Cooper pairs that form singlet states. Entangled 
states of electrons can thus be manipulated in transport experiments by extracting Cooper pairs out of a supercon- 
ductor. Several experiments using a superconductor as a source of entangled states of electrons have been proposed 
recently For instance it was shown in Ref. [1] that entangled states of electrons can be manipulated in a double dot 
experiment. A quantum teleportation experiment using three quantum dots has been proposed recently [2]. Another 
possible experiment has also been proposed in Ref. [3] in which a "beam splitter" is connected to a superconductor. In 
' this situation noise correlations can reveal information about electronic entanglement [3] . Other possible experiments 
\Q , in which several ferromagnetic electrodes are connected to a superconductor have been investigated theoretically in 
Refs. [4-8]. 

There is a rich physics occurring at a single ferromagnet / superconductor (FS) interface. For instance Andreev 
reflection is suppressed if the spin polarization of the ferromagnetic metal increases. This is because the incoming 
electron and the reflected hole belong to different spin bands. As a consequence Andreev reflection can occur only 
in the channels having both a spin-up and a spin-down Fermi surface [9]. This theoretical prediction was well 
verified in experiments [10,11] and it was shown that with high transparency interfaces the suppression of Andreev 
reflection by spin polarization can be used to probe the Fermi surface polarization [10]. The results of the Andreev 
reflection experiments compare well with another method based on spin polarized tunneling in the presence of Zeeman 
splitting [12]. Another phenomenon taking place at FS interfaces is that the pair amplitude induced in a ferromagnetic 
£j , metal can oscillate in space. This gives the possibility of fabricating tt junctions in which the Josephson relation 
changes sign [13-19]. It is also well established that FS multilayers present oscillations of the superconducting critical 
temperature as the thickness of the ferromagnetic layer is varied [20-25]. Other new phenomena related to the 
proximity effect have also been investigated in diffusive FS heterostructures [26-31]. 



Multiterminal hybrid structures consist of systems in which several spin polarized electrodes are connected to 
a superconductor and are controlled by crossed Andreev reflection processes in which the spin-up and spin-down 
electrons making the Cooper pair can tunnel in different ferromagnetic electrodes. Several theoretical predictions 
have been made. For instance the current circulating in one of the electrodes can be controlled by the voltage applied 
on another electrode [5,7]. The conductance can be described in terms of a conductance matrix [5] that can be 
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calculated from Keldysh formalism [7]. Other predictions concern the proximity effect in FSF trilayers [32]. It was 
shown in Ref. [8] that the two ferromagnetic electrodes of the FSF trilayer are coupled by pair correlations and 
that the superconducting order parameter in the ferromagnetic alignment can be larger than the superconducting 
order parameter in the antiferromagnetic alignment (Ap > Aaf)- O nc goal of our article is to find out the range of 
validity of this result and to bridge the gap with other predictions based on Usadel equations [32]. It was shown in 
the context of Usadel equations that the critical temperature (and therefore the superconducting order parameter) 
is larger in the antiferromagnetic alignment (Aaf > Ap) [32]. We show that the ballistic models that wc consider 
can lead also to Aaf > Ap if one takes into account the existence of a finite exchange field in the superconductor, 
in the spirit of Ref. [33]. There are thus two competing mechanisms for the proximity effect in FSF trilayers. The 
mechanism based on the exchange field involves the "11" and "22" diagonal elements of the Green's function of the 
superconductor connected to the ferromagnetic electrodes. The other mechanism based on pair correlations involves 
the "12" extra diagonal element. To understand in detail the mechanism based on pair correlations we investigate 
a device in which three ferromagnetic electrodes are coupled to a superconductor. As a simplifying assumption we 
suppose that the three electrodes are coupled to the same site in the superconductor. The fabrication of three contacts 
within a coherence length is not possible with present time state of the art technology. Nevertheless these systems 
with three electrodes provide idealized situations that are useful for understanding theoretically how the value of the 
superconducting order parameter depends on pair correlations induced in the ferromagnetic electrodes. 

From the point of view of the method we use two complementary approaches: (i) an analytical evaluation of the 
high energy behavior of the Gorkov function; (ii) exact diagonalizations of the Bogoliubov-dc Gcnnes equations. The 
analytical approach is based on the simplest ballistic models and can be a useful guideline for understanding more 
realistic systems in the future, such as multichannel systems or diffusive systems. An important difference between the 
analytical and numerical approaches lies in the fact that in the analytical calculation we suppose that the band-width 
of the superconductor is much larger than the superconducting gap like in realistic systems for which A/D ~ 10 -5 . 
As a consequence of this assumption the integral appearing in the self-consistency relation is dominated by the high 
frequencies. By contrast in the numerical simulation the band-width of the superconductor is not small compared 
to the superconducting gap (typically the ratio between the bandwidth and the gap is A/D ~ 1/5). The analytical 
calculation and the numerical simulation correspond to different regimes and this is why all the analytical results 
cannot be confirmed by the numerical simulation. The numerical simulation is useful to discuss the behavior of the 
local density of states (LDOS), the case of a partial spin polarization and the spatial variation of the pair amplitude. 

The article is organized as follows. Necessary technical preliminaries are given in section II. The mechanism based 
on the exchange field induced in the superconductor is investigated in sections III and IV. The remaining of the 
article is devoted to the other mechanism based on pair correlations. In section V wc discuss the proximity effect in 
FSF heterostructures and generalize the results obtained in Ref. [8] . Namely we show that the superconducting order 
parameter is larger if the ferromagnetic electrodes have a parallel spin orientation. Multiterminal structures in which 
three electrodes are connected to a superconductor are discussed in section VI. Friedel oscillations in the Gorkov 
function are analyzed in section VII. Concluding remarks are given in section VIII. 

II. PRELIMINARIES 
A. Green's function formalism 

In this section we provide the technical details on Green's function formalism that is used throughout the article. 

1. The models 

The superconductor is described by a BCS lattice model 

where A a is the pairing interaction. The brackets indicate that hopping is between nearest neighbors. The ferromag- 
netic electrodes are described by a lattice Stoner model 

Hstoncr = -t + ~ h J2 ( C i,1 Ci ^ ~ C ti°^) ' 
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where h is the exchange field. An exchange field smaller than the bandwidth corresponds to a partially polarized 
ferromagnet and an exchange field larger than the bandwidth corresponds to a half-metal ferromagnet in which only 
majority spins are present. We note pjf the normal state density of states of the superconductor and p% the spin-tr 
density of states in the ferromagnetic electrodes. We use Greek symbols a, (3, 7 for the sites in the superconductor and 
Latin symbols a, b, c for the sites in the ferromagnetic electrodes. The tunnel Hamiltonian coupling the superconductor 
and the ferromagnetic electrodes takes the form 

W = J^ia fc ,Q fc ( C a fc ,o- C ak,c + C a fc ,<r C afc ,"") > (1) 
k,cr 

where the sum over k runs over all contacts between the superconductor and the ferromagnetic electrodes. 



2. Green's functions 



The Green's functions of a connected system are obtained by solving the Dyson equation in the Nambu represen- 
tation: 



G R ' A — g R ' A + g R ' A (g) S ® G R ' A , 



(2) 



where the self-energy £ contains all couplings of the tunnel Hamiltonian given by (1). The Green's functions g 
correspond to the "disconnected" system in which t akiCtk = (sec Eq. 1). The Dyson equation (2) is used to calculate 
the Green's functions of the connected system in which an electron in the superconductor can make excursions in the 
ferromagnetic electrodes. 

We use the following notation for the Nambu representation of the advanced and retarded propagators of the 
disconnected system: 



f A > R (t,t<) g A - R (t,t>) 



with 



g A (t,t') = -ie(t-t')({c iA (t),c+ A (t')}) 

f A (t,t') = -ie(t-t')({c iA (t),c jd (t')}), 

and we use the following notation for the Nambu representation of the density of states: 



(3) 
(4) 



Pg Pi 
Pi Pg 

with p g = ^lm(g A ) and pf = ^Im(f A ). The Kcldysh Green's function is obtained through the Dyson-Keldysh 
equation 



I + G 1 



/ + E ® G A 



(5) 



where gf', = 2ii:n F (u> - Hi,j)pi,j ■ 



3. Superconducting order parameter and transport properties 

In equilibrium the Kcldysh Green's function defined by Eq. (5) simplifies into 

d+T - n F (uj - Mo ) (G A - G R ) , (6) 

where /io is the chemical potential, identical in all conductors. The form (6) of the Keldysh Green's function at 
equilibrium can be used to obtain the superconducting order parameter in the superconductor via the self-consistency 
equation 

A, = -^/^V'V). (7) 
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4- Green's function of an isolated normal metal 



Evaluating the spectral representation for a normal metal leads to 

A (,,\_ ma l "0 - U _ D A UR a jj 

VF 



fr 2 2nRa,fi 



exp[ —iUfR — i- 



The local density of states is given by 



pioc(w) = -Im [g A (u) 



1 man / , 1 



2tt 2 h 2 



vf 



(8) 



(9) 



5. Green's function of an isolated insulator 

We consider a semi-conductor model having a valence band at energy e(k) = —A — h 2 k 2 /(2m*) and a conduction 
band at energy e(fc) = A + h 2 h 2 /(2m*). This model behaves like a band insulator if \u\ < A. The Green's functions 
in the energy range \u>\ < A are given by 



9a A u ) = 



m*a 2 a 
h 2 2itR a . (: 



exp 



(-^^VATTJi?) - exp {-^^^~^ 



(10) 



6. Green's functions of an isolated superconductor 



The Green's function of an isolated superconductor with a uniform superconducting order parameter Ao is obtained 
by evaluating the spectral representation [7] : 



„ 2ma\ a 

MM = -fir- ^? exp 



R, 



a,f3 



{VA^ 



sin <p a ,ji 



-u) A 
A -co 



+ cos tp ai p 



-1 
1 



Above the superconducting gap the retarded Green's function takes the form 



9 a ,fi = h 2 2nR 



exp 



-u) A 
A —u 



+ COSip a .p 



-1 
1 



(11) 



(12) 



The phase variable in Eqs. (11) and (12) is <p ai p = kpR a ,i3- The "local" propagators corresponding to a = (3 are 
described by ip a ^ = fc^-Ro = n/2 [40]: 



5ioc(w) 



2ma 2 ao 



& 2irR ^A 2 - uj 2 



-oj A 
A -to 



(13) 
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FIG. 1. Schematic representation of the model considered in section II A 7 in which three half-metal ferromagnetic electrodes 
are connected to a superconductor. 
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7. Dyson matrix 



In this section we provide a derivation of the Dyson matrix in the simple case of three half-metal ferromagnetic 
electrodes connected to a superconductor (see Fig. 1). If A is an arbitrary site in the superconductor, the Dyson 
equation (2) becomes 



Qa n ,\ 



K l,l K l,2 



-K. 



2.1 



-K. 



2,2 



E 



K l,l ~ K 1,2 



f C 1 



a m ,X 



(14) 



with K^'° m - " a "' Q " 



= 9d *o B ,a„ffiJ' , and 



GQti-A /^id , A 
1,1 "1,2 

"2.1 "2,2 



(15) 



Eqs. (14) and (15) are valid for an arbitrary spin polarization. In the case of the heterostructure on Fig. 1 the explicit 
form of the Dyson matrix is the following: 



1 - K a a t 

1 "l l '-a, a 



-K h ' a t 



a, a 
c.a, 
1 1 £ a ,a 



t/3,b 



1,1 ''TjC 




"1,1 




1,1 


- K b,1 t 

1,1 7) c 




"1,1 




R b,\ 
1,1 






. G 14 . 




M 



In sections IV and V and in the Appendices we invert a similar form of the Dyson matrix in models involving four 
channels. 



B. Exact diagonalizations of the Bogoliubov-de Gennes Hamiltonian 

In the numerical simulations based on exact diagonalizations of the Bogoliubov-de Gennes Hamiltonian we use a 
lattice Hubbard model 

(a,0),cr a,a ct,a a.a a 

In Eq. (16) n QiCr = c+ a c a . a is the electron number operator at site a, fx is the chemical potential, h a cT = —ha z is 
the exchange field in the ferromagnetic region and a z — ±1 is the eigenvalue of the z component of the Pauli matrix. 
Vo is the on-site interaction. We use negative values of Vq corresponding to attractive interaction. To simulate the 
effect of depletion of the carrier density at the surface the site-dependent impurity potential is set to a sufficiently 
large value at the surface sites. This prohibits electron tunneling over these sites. Within a mean field approximation 
Eq. (16) reduces to the Bogoliubov-de Gennes equations [37,38]: 



£ A \ / u nA (r a ) \ _ f u nA (r a ) 



and 



such that 



£ A \ / u nd (r a ) \ _ ( u nA (r a ) 
A* -£) \v ni1 (r a ) J - en ^\v nA {r a ) 



(18) 



£u nt<T (r a ) = -f^y(r a +S)+ (/i 7 (r Q ) + /z) u„, CT (r a ) + h a a z u n ^{r a ) (19) 

Au n ^(r a ) = A (r a )u„ ;(T (r a ), (20) 
and where the pair potential is defined by 

A (r Q ) = Vb(c T (r Q ) Ci (r Q )). 
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In Eq. (19) 5 = x, —x, y, ~y denotes the four directions of the square lattice. 6 = z, —z corresponds to the two 
directions in the ferromagnetic electrodes. The self-consistency equation takes the form 



Ao(r a ) = = %1 



E 



u„^{r a )v* d (r a ) 



tanh ^ £ ^' 71 ^ +u„,|(r a )i;* ir (r a )tanh ^ P e ™rn 



(21) 



where (3 is the inverse temperature. We start from an approximate solution for the gap profile Ao(r Q ). After exact 
diagonalizations of Eqs. (17) and (18) we obtain u n ,a(r a ) and v n ^ a (r a ). The quasiparticle amplitudes are inserted 
into Eq. (21) and a new gap function is evaluated which is then inserted into Eq. (20) and we iterate until a sufficient 
precision has been obtained. Although the pair potential Ao(r a ) is zero in the ferromagnet the pair amplitude F{r a ) 
is non zero. The LDOS at site a is given by 



Pa {E) = -J2 [\u n Ar a )\ 2 .f'(E - e n ) + \v n , a (r a )\ 2 f'(E + e„)] , 

n, a 

where /' is the derivative of the Fermi function /(e) = i +cxp ^y^T) . 



(22) 



III. FRIEDEL OSCILLATIONS AT A SINGLE NF INTERFACE 



We start by considering a single interface between a normal metal and a spin polarized ferromagnet. We calculate 
the exchange field at site (3 in the normal metal, induced by the proximity of the ferromagnet (see Fig. 2). To order 



Normal metal 



Ferromagnet 



3 



a a 



FIG. 2. Schematic representation of the single-channel junction between a normal metal and a ferromagnet. R is the distance 
between sites a and (3. 



t 2 the Green's function at site (3 is given by 



Gb,B = 3/3, p + gfi,ata,aga,ata,ag a ,fi- 



(23) 



9a 8 i s gi yen by Eq. (8) and g^ a is decomposed in a real and imaginary part: g^ a = gi R a + igal- If the exchange 
field is small compared to the typical energy scales (being the energy band-width for a metal or the charge gap for 
an insulator) then the spin-up and spin-down Green's functions take the form 



diO 



where gi°a(w) is the Green's function in the absence of an exchange field. We deduce the existence of Friedel oscillations 
in the spin polarization at the Fermi surface: 



~(T) _ ~U) 
Pp,a P/3,3 

~(T) , -U) 

Pf3,P + P/3,0 



2nti 



1 / 1 



v F \2irR 



da [R) da {1) 
sin (2k F R) -p^ (0) + cos (2k F R) (0) 



du> 



du> 



(24) 



ferromagnetic electrode 




FIG. 3. The geometry used in the numerical simulations in section IV. A one-dimensional ferromagnetic electrode is con- 
nected to a two-dimensional superconductor. 
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FIG. 4. (a) Pair amplitude for different values of the exchange field for several sites in the superconducting region at the 
interface of one-dimensional ferromagnetic electrode and a two-dimensional superconducting system. Sites a, 7, f3 belong to 
the superconductor while the sites ci , C2 belong to the ferromagnetic electrode, (b) Pair amplitude for several sites in the 
ferromagnetic region for different values of the exchange field. The pair amplitude in the superconductor does not vary much 
with the exchange field. The pair amplitude in the ferromagnetic electrode shows an oscillatory behavior. The pair amplitude 
decays monotonically in the absence of spin polarization in the ferromagnetic electrode. 



For an insulating ferromagnet we obtain from Eq. (10) that dg a ,a /dui(0) is negative which leads to a positive spin 
polarization at site a. For a metallic ferromagnet we obtain from Eq. (9) that dgi)a /^(0) is positive. Since we 
use kpRa = Tt/2 we deduce that spin polarization at site a is negative. As far as the magnitude of the effec- 
tive exchange field /i c ff is concerned we see that h c s is of order (t a ^ a /ep) 2 h cx in the metallic case and of order 
t 2 , Q A -1 / 2 (m*aQ//i 2 ) 3 ' /2 /i ex in the insulating model. 

We note that in the insulating case the limit R — > is well-defined and the value of spin polarization at site a does 
not depend on the value of the ultraviolet cut-off i?o introduced at short length scale. In the metallic model the value 
of the spin polarization is diverging if one takes the limit R — > 0. Therefore the sign of spin polarization at site a 
depends on the value of the short distance cut-off Rq- In the analytical approaches we use a regularization in which 
we introduce a cut-off equal to kpRo = 7r/2 that is chosen similarly to Ref. [40]. The justification of this choice of Rq 
is that it leads to a physically acceptable value of the local density of states. The value of the density of states at the 
Fermi energy for R — Rq is given by 

man an 

AT 2ttR 

With the choice k F R = n/2 of the short distance cut-off the local density of state is equal to p a , a {Ro) = 
ma^kp/ (ir 2 h 2 ), a positive value that is not far from the exact density of states obtained for Rq = (see Eq. (9)). 



IV. PROXIMITY EFFECT AT A SINGLE FS INTERFACE 



A. Local density of states 

We first use exact diagonalizations of the Bogoliubov-de Gennes equations (see section II B) to describe the proximity 
effect at a single FS interface. This constitutes a test of the numerical method and we recover some of the results for 
the LDOS obtained recently by means of a recursion method (see Ref. [36]). 

We consider a two-dimensional superconducting system of 30 x 30 sites and we suppose fixed boundary conditions 
by setting the impurity potential yJ = lOOt at the surface. The temperature is ksT — 0. It and the local attractive 
interaction in the superconducting region is Vq = — 3.5t. On top of the superconductor is attached a one-dimensional 
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FIG. 5. (a) The LDOS at site 7 in the superconductor for different values of the exchange field, for the geometry of Fig. 4. 
(b) The LDOS at site ci in the ferromagnetic electrode for different values of the exchange field. The subgap LDOS is due to 
Andrccv bound states. 



superconducting electrode of 50 sites (see Fig. 3). The transparency of the interface is controlled by changing the hop- 
ping element t c connecting sites on both sides of the interface and we restrict here to the case where the transparency 
of the interface is the same as inside the superconductor. In this case the pair amplitude in the superconducting 
system is not really modified as seen in Fig. 4(a) while in the ferromagnetic region the pair amplitude oscillates 
around zero and the period of oscillations decreases with increasing the exchange field (see Fig. 4(b)). In the case of 
a zero exchange field the pair amplitude is decaying monotonically in the normal metal. 

Due to the proximity effect the LDOS shows a gap structure even for the sites within the ferromagnet (see Fig. 5(b)). 
The conductance peaks within the gap are due to Andreev bound states [39] and have been discussed recently for a 
three dimensional FS interface using a recursion method [36]. The residual values of the LDOS are reduced by the 
increase of the exchange field. The Andreev bound states move towards the Fermi level and cross the Fermi level with 
increasing the exchange field. 



B. Exchange field in the superconductor 



From the discussion in section III we deduce that a finite exchange field is induced in the superconductor [33] . The 
effect was also found in recent numerical simulations in Ref. [41]. In the metallic case and in the tunnel limit the 
magnitude of the exchange field is of order (t a , a /^F) 2 h c ^. 

An exchange field in a superconductor is a pair-breaking perturbation that tends to dissociate Cooper pairs [37,42] . 
In this situation we expect that Cooper pairs couple to the exchange field averaged over a length scale equal to the 
BCS coherence length. A qualitative argument can be made from the NF model in section III. Since the discussion 
is similar in the two cases of metallic and insulating ferromagnets we consider a metallic ferromagnet model only. We 
obtain 

~(T) i \ ~U) t \ *«,<A X ( ma V\ ( a o\ 2 rQ /, . / \i 



from what we deduce 

r+oc r D 



/>-t-oo C r 

/ 4nR 2 dR / ^(w)-p^M 

JRn JO L 



iRo JO L " ' " J 37T 3 \ H 2 J j 3kpRo 

With kpRo — tt/2 the spin polarization induced at long distance in the superconductor is positive but the local spin 
polarization at the contact is negative. 
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The reduction of the superconducting order parameter due to pair-breaking effects was observed in FSF trilayers 
with insulating ferromagnets in Refs. [43,44]. The effect is not visible in the point contact experiments in Ref. [10] even 
though one might expect a reduction of the superconducting gap at the interface that would increase with increasing 
the exchange field. It is likely that the point contact geometry used in Ref. [10] plays a role and that the effect should 
be more pronounced for extended contacts. 




(a) (b) 

FIG. 6. The two models considered in sections V B and VII. In (a) the two ferromagnetic electrodes are connected to the same 
site in the superconductor. In (b) the two ferromagnetic electrodes are connected to two different sites in the superconductor. 
The two ferromagnetic electrodes contain a spin-up and a spin-down conduction band. The models are thus solved in a 4 x 4 
formalism. 
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FIG. 7. The model considered in section VI. In (a) the three electrodes are connected to the same site in the superconductor. 
In (b) the three electrodes are connected to two different sites in the superconductor. We use the model (a) to calculate the 
superconducting order parameter . 



V. PROXIMITY EFFECT IN FSF HETEROSTRUCTURES: LOCAL MODELS 



Now we reconsider the proximity effect in a FSF heterostructure in which two ferromagnetic electrodes are connected 
to a superconductor. In this section as well as in section VI we consider "local" models in which two or three 
ferromagnetic electrodes are connected to the same site in the superconductor (see Figs. 6-(a) and 7-(a)). 



A. Motivation 



Two effects can play a role in the determination of the self-consistent order parameter: 
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FIG. 8. Variation of F(t,T) for different values of T < 1/2. We obtain F(t,T) > which shows that the ferromagnetic 
superconducting order parameter is larger than the antiferromagnetic superconducting order parameter for large interface 
transparencies. The microscopic interaction U is such that the bulk superconducting gap is A bu ik = 1 and the bandwidth of 
the superconductor is D = 100. 

(i) The proximity effect that takes place already at a single NS interface. A pair amplitude is generated in the 
normal metal and there is a reduction of the order parameter on the superconducting side of the interface [37,42]. 
For the models on Fig. 6 the two electrodes ending at sites "a" and "b" are coupled by pair correlations only in 
the antiferromagnetic alignment. The case of an antiparallel spin orientation is thus qualitatively similar to a 
NS interface in which case the order parameter in the superconductor is reduced at the interface. On the other 
no pair correlations are generated among the ferromagnetic electrodes in the parallel alignment if we consider 
half-metal ferromagnets. The case of parallel magnetization is thus similar to an isolated superconductor, from 
what we deduce that the superconducting order parameter is larger in the ferromagnetic alignment. 

(ii) Pair breaking effects due the finite exchange field in the superconductor. The exchange field is a pair-breaking 
perturbation that tends to reduce the superconducting order parameter. 

A relevant question is to determine which effect would control real experiments on metallic FSF trilayers. It is 
likely that the answer depends on the geometry of the devices. For small interface transparencies the exchange 
field in the superconductor is of order (t/e F ) 2 (see section IV). For a thin film in a parallel field the pair breaking 
parameter is equal to a = De 2 H 2 d 2 / (Qhc 2 ) with D the diffusion constant, H the magnetic field in Tesla, d the width 
of the superconductor. The difference between the critical temperatures in the ferromagnetic and antiferromagnetic 
alignments is given by fcsAT c = 7ra/4 [42] which is proportional to t 4 because H is proportional to t 2 . On the other 
hand if the physics is dominated by the proximity effect we deduce from Rcf. [6] that the difference between the 
superconducting order parameters in the ferromagnetic and antiferromagnetic alignments is of order Af — Aaf = 
t 4 A(°) /(Ue%), with A<°) = Dexp[-l/(Up N )] the BCS superconducting order parameter. We see that AT C due to 
pair-breaking effects is proportional to d 2 and we deduce that pair-breaking effects are reduced if the width of the 
superconducting layer is reduced. 

In the following we focus on the effect of pair correlations induced in the ferromagnetic electrodes and we take for 
granted that the width of the superconductor is small enough so that pair-breaking effects can be neglected. We make 
the further simplifying assumption of considering that all ferromagnetic electrodes are connected to the same site and 
we calculate the superconducting order parameter at this site. The case of two electrodes at a finite distance will be 
considered in section VII. The qualitative physics occurring in the schematic model with all electrodes connected to 
the same site can be described the following rule: increasing pair correlations in the ferromagnetic electrodes tends 
to reduce the value of the self- consistent order parameter in the superconductor. 

To illustrate this let us first consider the case of two ferromagnetic electrodes connected to a superconductor (see 
Fig. 6- (a)). Pair correlations among the two ferromagnets are stronger in the antiferromagnetic alignment than in the 
ferromagnetic alignment. As a consequence of the rule given above the superconducting order parameter is larger in 
the ferromagnetic alignment. To make a more refined test of this rule we consider in section VI a system in which three 
electrodes are connected to a superconductor (see Fig. 7- (a)). For simplicity we consider that the three ferromagnetic 
electrodes are half-metal ferromagnets. Electrode c is supposed to have a spin-down orientation. Electrode a has a 
spin-up orientation. Electrode b can have a spin-up orientation (ferromagnetic alignment) or a spin-down orientation 
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(antifcrromagnetic alignment). If we consider first that t Cj7 is small we see that pair correlations are formed mainly 
among electrodes a and b if electrodes a and b have an opposite spin orientation. We deduce that Ap > Aaf if t Cil 
is small. We consider now that t cn is large and electrode c has a spin-down orientation. If electrodes a and b have a 
spin-up orientation we see that pair correlations are formed mainly among electrode a and c and among electrodes b 
and c. If electrode a has a spin-up orientation and electrode b has a spin-down orientation we see that pair correlations 
are formed mainly among electrodes a and c. As a consequence the proximity effect is stronger if electrodes a and b 
have a parallel orientation. We deduce from the rule given above that Ap < Aaf if t Cil is large. 



electrodes 




Y 

FIG. 9. The geometry used in the numerical simulations of the Bogoliubov-de Gennes Hamiltonian. Two one-dimensional 
ferromagnetic electrodes are connected to the same site of a two-dimensional superconductor. Electrode ending at site "a" is 
represented by filled circles and electrode ending at site "b" is represented by filled squares. 




FIG. 10. Difference between the LDOS in the ferromagnetic and antiferromagnetic spin alignments for the geometry on 
Fig. 9. The LDOS is represented as a function of E/t, for two values of the exchange field (h = 1 and h — 2). The upper panel 
corresponds to the LDOS at site 7 and the lower panel corresponds to the LDOS at site a (see Fig. 9). The subgap LDOS is 
larger in the ferromagnetic alignment. 
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FIG. 11. The difference between the pair amplitudes in the ferromagnetic and antiferromagnetic alignments for two electrodes 
connected to the same site 7. 

B. Two ferromagnetic electrodes connected to the same site 

1. Sign of Af — Aaf 

Let us start with the situation where two ferromagnetic electrodes are connected to the same site (see Fig. 6- (a)). 
We suppose that each ferromagnetic electrode is partially spin polarized and contains two spin channels. The model 
is thus solved in a 4 x 4 formalism. The Green's function G^'a 2 is found to be 

Ga'a 2 = -™p$- - ™p0-^ {-/Kt) - f( x a,l) ~ /Kt) " /Kl) + V{x aA )f{x bA ) + 2f(x ad )f(x bd ) (25) 
U) CO D 

+ f(Xa,l)f(Xa,l) + f( x a,f)f( X b,l) + f( x a,l)f( x b,l) + f( x b,f)f( x b,l) ~ f( x a,f)f( X a,l)f( X b,l) 

~ f(Xa,l)f(Xa,l)f( x b,l) - f( x a,l)f( x bA)f( x b,l) ~ f( x a,l)f( x b,-[)f( x b,l) - f( x a,-[)f( x a,l)f( x b,-[)f( x b,l)} > 

with 

V = (1 - f(x ad )f(x bd )) (1 - /(x a , T )/(^, T )) . (26) 

The interface transparencies are parametrized by f{x{) — Xi/(1 + x{), with x ad — ^ 2 t^ a p^ j x a,i — n2 ta, a Pa 1P0 > 
£b,T = n2 t 2 pPb ]Po j x b,l — ^t 2 pPb iPo ■ The values of f{x{) are such that < f(xi) < 1/2. We do not consider the 
regime 1/2 < /(a;,) < 1 for which we find that the self-consistency relation is unstable. 

We suppose that the two ferromagnets have an identical density of states and that the two contacts have identical 
transparencies. In the ferromagnetic alignment we use the notation f(x ad ) — f{x bd ) — T and f{x ad ) = f(x bd ) = t. 
In the antiferromagnetic alignment we use the notation f(x ad ) — f{x bd ) — T and f{x ad ) = f{x bd ) = t. To second 
order in T and t we obtain 

G^f-G^ F = -i^(T-tf. (27) 

As a consequence in the regime of small interfaces transparencies we obtain Ap > Aaf- The regime of a large interface 
transparencies can be treated by evaluating numerically the difference between the Gorkov functions „ F — G]; 2 ^ AF = 
-iirp$jjF(T,t). It is visible on Fig. 8 that F(T,t) is positive for any value of T < 1/2 and t < T so that the 
ferromagnetic superconducting order parameter is larger than the antiferromagnetic superconducting order parameter, 
in agreement with Ref. [8]. 

C. Local density of states and pair amplitude 

We calculated the LDOS for a partial spin polarization by means of exact diagonalizations of the Bogoliubov-de 
Gennes Hamiltonian. The geometry of the simulation is represented on Fig. 9. It is visible on Fig. 10 that the low 
energy LDOS is larger in the ferromagnetic alignment than in the antiferromagnetic alignment. 
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The variation of the difference Fp — Faf between the pair amplitudes in the ferromagnetic and antiferromagnetic 
alignments as a function of the exchange field is shown on Fig. 11. The numerical simulation coincides with the 
analytical model for a large exchange field in the sense that the pair amplitude is larger in the ferromagnetic alignment. 
However for small values of the exchange field we obtain Fp > Faf or Fp < Faf depending on which lattice site is 
considered. The possible origin of the discrepancies between the numerical simulation and the analytical model are 
analyzed in the concluding section. 

VI. PROXIMITY EFFECT IN MULTITERMINAL HYBRID STRUCTURES: LOCAL MODELS 

In this section we consider an heterostructure in which three electrodes are connected at the same site to a supercon- 
ductor. Electrodes a and b on Fig. 7 will be considered to be half-metal ferromagnets so that the models are solved in 
a 4 x 4 formalism. Using partially polarized ferromagnets would require to solve a 6 x 6 formalism for which we could 
not carry out the analytical calculation in the regime of large interface transparencies. The case of partially polarized 
ferromagnets will be treated numerically within exact diagonalizations of the Bogoliubov de Gennes equations. 

A. Three electrodes connected at the same site 

We consider a model in which three electrodes are connected to the same site a (see Fig. 7-(a)). We determine the 
variation of Ap — Aaf as a function of the transparency t c of the contact with electrode c. We show that there is no 
change in the sign of Ap — Aaf as t c increases if electrode c is a normal metal. There is a change of sign in Ap — Aaf 
as t c increases if electrode c is a ferromagnet with a spin-down orientation. This is in agreement with the qualitative 
discussion in section VA. 

B. Sign of A F - Aaf 

In the ferromagnetic alignment the inversion of the 4x4 Dyson matrix leads to 

Gi' 2 Q = ~i*p S Q - - i^-^r {-/(*„) - f(x b ) - f(x cA ) - f(x cd ) + 2f(x a )f(x b ) + 2f(x a )f(x cA ) (28) 
+ 2f(x b )f(x Ci i) + f(x a )f(x cd ) + f(x b )f(x cd ) + f(x cA )f(x cd ) - 3f(x a )f(x b )f(x cA ) - f(x a )f(x b )f(x cd ) 

- f(Xa)f(Xc,l)f(x c ,l) - f(Xb)f(x c ,l)f(x c ,l) + f( x a).f(x b )f(XcA).f(Xc,l)} > 

with 

V ¥ = 1 - f(x a )f(x b ) - f{x a )f{x cA ) - f{x b )f{x cA ) + 2f(x a )f{x b )f(x cA ). 
In the antiferromagnetic alignment we find 

G^l = -inp$- - ittpI--^- {-f(x a ) - f(x b ) - /(x c , T ) - f(x Ctl ) + 2f(x a )f(x cA ) + 2f(x b )f(x c>l ) (29) 

+ f(x a )f(x b ) + f(x a )f(x cd ) + f(x b )f(x cA ) + f{x cA )f(x cd ) - f(x a )f(x b )f(x cd ) 

- f(x a )f(x b )f(x cd ) - f(x b )f(x cA )f(x cd ) - f(x a )f(x cA )f(x cd ) - f(x a )f(x b )f(x c ,t)f(x c ,i)} , 

with 

X>AF = 1 - f(x b )f(x Ctl ) - f{x a )f(x cd ) + f(x a )f(x b )f(x cd )f(x Ctl ). 

The variation of A£ a - A£ F a as a function of the transparency of the contact with electrode c is shown on Fig. 12. 
It is visible on this figure that Ap > Aaf for all values of the transparency f(x c ) with interface c. 

We also evaluated Ap — Aaf in the case where electrode c is a half-metal ferromagnet with a spin-down orientation. 
We obtain a change of sign in A^ a — as the transparency of the interface with the ferromagnetic electrode c is 
increased (see Fig. 13). This change of sign can be described by expanding the Gorkov functions to second order in 
f(x a ) and f(xb): 
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FIG. 12. Variation of A^ a — as a function of the transparency of the contact with the normal electrode f(x c ^) = f(x Ci i) 
for different values of the transparency of the contacts with the ferromagnets f(x a ) = f(xb)- The different curves correspond to 
(a) f( Xa ) = f{x b ) = 0.05, (b) f( Xa ) = f(x b ) = 0.1, (c) f( Xa ) = f(x b ) = 0.15, (d) f( Xa ) = f(x b ) = 0.2, (e) f(x a ) = f{x b ) = 0.25. 
We obtain Af > Aaf- The microscopic interaction U is such that the superconducting order parameter of the isolated 
superconductor is Abuik = 1 and the bandwidth of the superconductor is D = 100. 



£tl,2,Ferro _ „•„ „S ^ „-_„S^ 



G. 



l,2,Antiferro 



™Po "TPo — [~f(x a ) - f(x b ) - f{x c ,i) ~ f(x c ,i) + 2 f(x a )f{x b ) + 2/(x a )/(x c , T ) 

UJ UJ 

+ 2f{x b )f{x cA ) + f(x a )f(x cd ) + f(x b )f(x cd ) + f(x cA )f(x cd )} 

-mp%- - inp$- [-f(x a ) - f{x b ) - f(x cA ) - f{x cd ) + 2f(x a )f{x cA ) + 2f(x b )f(x cd ) 

+ f(Xa)f(x b ) + f(x a )f(x cd ) + f{x b )f{x cA ) + f(x c ,l)f(x cd )] , 



(30) 



(31) 



from what we deduce that the difference between the ferromagnetic and antiferromagnetic superconducting order 
parameters changes sign if f(x C) i) — f(x a ) + f(x c ^), in agreement with Fig. 13. 



C. Three node circuit model 



The presence or absence of sign changes in — A^ discussed in section VI B can also be obtained in the 

framework of a phenomenological circuit model [6] . In this model the superconducting order parameter is given by 



A = D exp 



Up 



N 



(l + 7rp A rr T )(l + 7 rp A rr i ) 



where pn is the density of states in the normal state, and where T a is the total spectral line-width of spin-cr electrons, 
obtained as the sum of the spectral line- widths associated to each electrode [6] : 

k 

where ts, ak is the tunnel amplitude connecting the superconductor and the ferromagnetic electrode a/.. 

Electrodes a and b are supposed to have the same spin polarization. We denote by 7 = |£s ja | 2 Pa,T the spectral line- 
width associated to majority-spin electrons in electrode a and by A7 = \ts. a \ 2 Pa,i the spectral line-width associated 
to minority-spin electrons in electrode a. The spin polarization of the ferromagnetic electrodes is P = (1 — A)/(l + A). 

We first suppose that electrode c is a ferromagnet having a spin-down orientation with a spectral line-width 70 
associated to majority-spin electrons and a spectral line-width A70 associated to minority-spin electrons. If the 
ferromagnets a and b have a parallel spin orientation the total spectral line-widths are given by = 2-f + A70 and 
T[ = 2A7 + 70. If the ferromagnets a and b have an antiparallel spin orientation the total spectral line- widths are 
given by = (1 + A)7 + A70 and T| = (1 + A)7 + 70. We obtain 



^AF 



exp 



k 2 Pn 
U 



(l-A) 2 (7 2 -7o 2 ) 
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FIG. 13. Variation of A^ a — as a function of the transparency of the contact with the spin-down ferromagnetic electrode 
f(x c ,i) (f(xc,t) = 0) for different values of the transparency of the contacts with the ferromagnets f(x a ) = f(x b )- The different 
curves correspond to (a) f(x a ) = f(x b ) = 0.05, (b) f(x a ) = f(x b ) = 0.1, (c) f(x a ) = f(x b ) = 0.15, (d) f(x a ) = f(x b ) = 0.2, (e) 
f(xa) = f(x b ) — 0.25. Af — Aap changes sign as the transparency of the contact with electrode c increases. The microscopic 
interaction U is such that the superconducting order parameter of the isolated superconductor is Abuik = 1 and the bandwidth 
of the superconductor is D — 100. 
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from what we deduce that Ap > Aaf if 7o < 7 and Ap < Aaf if 7o > 7> i n agreement with the microscopic model 
discussed in section VI B. A similar calculation in the case where electrode c is a spin-up ferromagnet shows that 
there is no sign change in Ap — Aaf- 

The presence of a sign change in Ap — Aaf as 7 increases was anticipated in section VA. The value 7 = 70 at 
which the change of sign takes place can be understood from the following symmetry argument. In the ferromagnetic 
alignment of electrodes a and b, electrodes a and b have a spin-up orientation and electrode c has a spin-down 
orientation. In the antifcrromagnctic alignment electrode a has a spin-up orientation and electrodes b and c have a 
spin-down orientation. Therefore if the three contacts have the same transparency (namely if 7 = 70) we see that the 
case of a ferromagnetic alignment of electrode a and b can be deduced from the case of an antifcrromagnctic alignment 
just by reversing the spin of the three ferromagnetic electrodes. As a consequence we obtain Ap — Aaf = if 7 = 70 
and therefore a sign change in Ap — Aaf if 7 = 7o- This argument is also in agreement with the fact that Ap — Aaf 
changes sign for f(x c ) = f(x a ) = f(xb) in the microscopic model in which three electrodes are connected to the same 
site (see Fig. 13). 

We suppose now that electrode c is a normal metal with a spectral line- width 70 associated to spin-up and spin-down 
electrons. We obtain 



A F 

Aa7 = cxp 



which leads to Ap > Aaf for all values of 70, in agreement with section VI B. 



D. Diagonalizations of the Bogoliubov-de Gennes Hamiltonian 

We consider exact diagonalizations of the Bogoliubov-de Gennes Hamiltonian in a situation where a two-dimensional 
superconducting system is connected to two one-dimensional ferromagnetic electrodes and to a one-dimensional normal 
metal electrode (see Fig. 14). The proximity effect in the normal metal electrode can be controlled by the spin 
orientation of the two ferromagnetic electrodes in the sense that the LDOS in the normal metal electrode depends on 
the relative spin orientation of the ferromagnetic electrodes. 

The LDOS at site c is shown on Figs. 15(a) and 15(b) in the cases of antiparallel and parallel spin orientation of 
electrodes a and b. It is visible that for a sufficiently large exchange field the LDOS in the normal metal electrode 
depends on the relative spin orientation of the ferromagnetic electrodes. As a consequence it may be possible to 
control the proximity effect in the normal metal by changing the spin orientation of the ferromagnetic electrodes. 
This may be of interest in view of STM experiments similar to Ref. [46]. 
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FIG. 14. The geometry used in the numerical simulation in which three one- dimensional electrodes are connected to a 
two-dimensional superconductor. Electrode ending at site "a" is ferromagnetic and is represented by filled circles. Electrode 
ending at site "b" is ferromagnetic and is represented by filled squares. Electrode ending at site "c" is a normal metal and is 
represented by filled triangles. 



VII. FRIEDEL OSCILLATIONS IN THE GORKOV FUNCTION 



Sections V and VI were devoted to the analysis of "local" models in which several electrodes are connected to 
the same site in the superconductor. Now we consider the same problem but with a finite distance between the 
ferromagnetic electrodes. As discussed in section VA we suppose that the width of the superconductor is small 
enough so that pair-breaking effects due to the finite exchange field can be disregarded. We consider the geometry on 
Fig. 6 in which a superconductor is in contact with two ferromagnetic electrodes a and a'. Compared to Fig. 6 we use 
the notation a' and a' rather than (3 and b for the second electrode. We keep the notation f3 for an arbitrary site in 
the superconductor at which we calculate the superconducting order parameter . We suppose that the two electrodes 
a and a' are at a distance much larger than the Fermi wave-length. 



A. Form of the Gorkov function 

Electrodes a and a' are characterized by the microscopic propagators g a , a and g a ',a'- To order t 4 the Green's 
function at site (3 is given by 

G/3,/3 = g/3,/3 + gi3,ata, a ga,ata,aga,l3 + g/3,a'ta' ,a> <?a' ,a't a ' ,a> (la' ,j3 (32) 
~t~ g@,ata,aga,ata.aga.ata.aga.ata,aga,j3 

(33) 

gft,ata.aga.ata.aga.a'ta' ,a' §a f ,a'ta' ,a f jja' (34) 
~t~ g/3,a'ta' .a' .a'ta' ,a' §a f ,ata,aga,a^a,aga,j3 (35) 
+ gf3,a'ta',a'ga',a'ta',a'ga',a'ta',a'ga',a'ta',a'ga',l3- (36) 

To lowest order the information about the exchange field in the superconductor is contained in the "11" and "22" 
components of the two terms in Eq. (32). The "12" component of the two terms in Eq. (32) contains the information 
about the effect of a single interface on the Gorkov function at order t 2 . The two terms in Eqs. (33) and (36) describe 
a single interface to order t 4 and are not directly relevant to our discussion since they are just a small correction of the 
terms in Eq. (32) obtained at order t 2 . The "12" component of the two terms (34) and (35) describe the modification 
of the Gorkov function due to processes coupling the two interfaces. Each of the two terms (34) and (35) gives rise to 
three terms in the Gorkov function when the summation over spin indexes is made explicit: 

( 34 ) = ^:^a,a5a;a^a,a5a;a' ^ "^a'3a' 2 a' ^ a^a'35 2 /3 ( 37 ) 

, -Ll f -1,1, -1,2 , -2,2 , -2,2 /oo\ 

' 9f3,ct a , a 9a,aJ'a,,aga,a' T ' a ', a '9a l ,a' T ' a ', a '9a',f3 \y°) 

, —1,1 i -l,lj- -1.1 j- -1,1 j- -1,2 fOQ\ 

9/3,a ta , a 9a,a t a,aga,a' t a' ,a' g a / , a ^a' ,a' g a \/3 ) 

(35) = .g^,^ Q ^a'3a' 2 a' ^ a^a'5 2 t ^ Q ^a,a3a;a ^ a,a5a, 2 3 (40) 
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FIG. 15. (a) The LDOS as a function of E/t at the interface of a multiterminal junction made of two one-dimensional 
ferromagnetic electrodes and a one-dimensional normal metal electrode that are connected on top of a two dimensional su- 
perconducting system. The LDOS is calculated at site c in the normal metal electrode. The upper panel corresponds to the 
antiferromagnetic alignment and the lower panel corresponds to the ferromagnetic alignment. 
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(41) 
(42) 



B. Friedel oscillations in the Gorkov function 

Now we use a single channel model to discuss Friedel oscillations in the Gorkov function. The terms of order t 2 are 
discussed in section VII Bl and the terms of order t 4 are discussed in section VII B 2. 

1 . Terms of order t 2 

Let us start with the terms of order t 2 . Since these terms do not couple the two interfaces we consider a single 
interface only. We obtain easily the expression of the Gorkov function to order t 2 : 



r i,2.A _ 2ma\ a A 
M ~ 1 H 2 2ttR Vc^ 2 - A 2 



( /,-' 



2Wuj 2 - A 2 R 



^u 2 - A 2 



(43) 
(44) 



If site (3 is equal to site a we use the regularization UfRq = 7r/2 discussed in section III to obtain the self-consistent 
order parameter at the interface: 



2D exp 



1 f 2771a 2 , a 
U \ nh 2 2vri? 



2ma?, a,Q 2 
1 H 2 2.R^ {P ^ +P ^ ) 
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FIG. 16. Variation of F(R) defined by Eq. (45). We used the parameters k F = 1, A = 1, D = 10 5 and v F = D/k F . The 
spin-up density of states is p Ql -f = 1. The spin-down density of states is p a ,i = (a), p a ,i = 0.5 (b) and p a .\ = 1 (c). 



As expected the superconducting order parameter at the interface is reduced because of the proximity effect. 
To discuss the proximity effect in the bulk of the superconductor we write the pair amplitude under the form 



f 

J A 



Im 



G 



1,2, A 

0,0 



dio = — 



L 



A 



2ma 2 ao 
K 2 2irR J A ^lo 2 - A 2 



du> + 



2ma 2 ) 



with 



F(R) 



sin [kpR] 

i?2 



/' 

J A 



Pa,] sin ■ 



k F + 



2V^ 



VF 



-R 



+ Pad sin • 



kF 



A 2 



VF 



-R 



^LO 2 - A 2 



duj. (45) 



The variations of F(R) as a function of R are shown on Fig. 16. It is visible that F(R) oscillates in space but is 
always positive which means that the superconducting order parameter is reduced in the bulk of the superconductor 
and shows Friedel oscillations. 
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R 

FIG. 17. Variation of G(R) defined by Eq. (47). We used the parameters k F = 1, A = 1, D = 10 5 and v F = D/k F . The 
spin-up density of states is p a ,t = 1. The spin-down density of states is p a ^ = 0.2. The distance between the two interfaces is 
R a ,a' = 100 and the function G(R) is calculated as a function of R — R a ,p = R a , a ' — R a ',p- 



2. Terms of order t 4 

Now we consider the terms of order t 4 that couple the two interfaces, namely the terms (34) and (35). Using 
Eqs. (37) - (42) we obtain 
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(46) 

(47) 
(48) 
(49) 
(50) 
(51) 
(52) 



The variations of G(R) are shown on Fig. 17. It is visible that G(R) is larger in the antiferromagnetic alignment 
and therefore in the bulk of the superconductor the superconducting order parameter is larger in the ferromagnetic 
alignment. This shows that the effect is not a specificity of the local models analyzed in sections V and VI in which 
all electrodes are connected to the same site. On the contrary the effect occurs also in the models on Figs. 6-(b) 
and 7-(b) in which the distance between the ferromagnetic electrodes is large compared to the Fermi wave-length. 



VIII. CONCLUSION 

To conclude we have provided a detailed investigation of the proximity effect in multi-connected hybrid structures 
in which several electrodes are connected to a superconductor. We have pointed out the existence of two mechanisms 
involved in the determination of the superconducting order parameter . One is related to the existence of an exchange 
field in the superconductor that was first pointed out in Ref. [33]. The other mechanism takes place already at 
a single NS interface. We have first reconsidered FSF heterostructures in which two ferromagnetic electrodes are 
connected to a superconductor at a distance smaller than the superconducting coherence length. We have found 
that within a "local" model the superconducting order parameter in the ferromagnetic alignment is larger than the 
superconducting order parameter in the antiferromagnetic alignment, in agreement with a different model discussed 
in Ref. [8]. We have shown that the zero-energy LDOS in the ferromagnetic alignment is larger than the zero-energy 
LDOS in the antiferromagnetic alignment. In the case where two ferromagnetic and a normal metal electrode are 
connected to a superconductor we have found that the LDOS in the normal metal depends on the spin orientation in 
the ferromagnetic electrodes. If the normal metal is replaced by a ferromagnetic metal with a spin-down orientation 
we have found that Aaf > Ap for high transparencies and we have provided two analytical models for this behavior 
(an inversion of the 4x4 Dyson matrix for half-metal ferromagnets and another approach based on a circuit model 
for partially polarized ferromagnets). We have pointed out that this behavior could be understood from a simple rule 
stating that for this model the increase of pair correlations in the ferromagnetic electrodes generates a reduction of 
the superconducting order parameter at the interface. This behavior is confirmed by the exact diagonalizations of the 
Bogoliubov-de Gennes equations for large values of the exchange field but not for small values of the exchange field. 

One explanation to the discrepancy between the two models is that in the Green's function approach only the 
contribution of energies far above the superconducting gap has been taken into account. This is legitimate if the 
band-width D of the superconductor is much larger than the superconducting gap (in real systems one has typically 
D/A ~ 10 5 ). This is justified by the fact that the integral in Eq. (7) is dominated by the high-energy behavior of the 
Gorkov function G^"'^"' 1,2 (w) ~ 1/lo. On the other hand in the exact diagonalizations of the Bogoliubov-de Gennes 
Hamiltonian we used small system sizes and replaced the right hand side of Eq. (7) by a sum over all energy levels. 
It is thus not surprising that the two procedures (namely the high-energy behavior for large systems and a sum over 
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all energy levels for small systems) can lead to different results in some cases. Another possible source of discrepancy 
between the two approaches lies in the fact that the ferromagnetic electrodes are one-dimensional and have a finite 
size in the exact diagonalizations whereas the ferromagnetic electrodes are described by the local density of states of 
a three dimensional metal in the Green's function approach. 



ACKNOWLEDGEMENT 



The authors acknowledge fruitful discussions with H. Courtois, D. Feinberg, M. Giroud. 



[1] M. S. Choi, C. Bruder and D. Loss, Phys. Rev. B 62, 13569 (2000); 

P. Recher, E. V. Sukhorukov and D. Loss, Phys. Rev. B 63, 165314 (2001). 
[2] D. Feinberg, T. Martin and O. Sauret, cond-mat/0203215. 
[3] G.B. Lesovik, T. Martin, and G. Blatter, arXiv:cond-mat/0009193. 
[4] G. Deutscher and D. Feinberg, App. Phys. Lett. 76, 487 (2000). 
[5] G. Falci, D. Feinberg, and F.W.J. Hekking, Europhysics Letters 54, 255 (2001). 
[6] R. Melin, J. Phys.: Condens. Matter 13, 6445 (2001); 

R. Melin, to appear in the Proceedings of the XXXVIth Rencontres de Moriond, T. Martin and G. Montambaux Eds., 

EDP Sciences (2001). 
[7] R. Melin and D. Feinberg, Eur. Phys. J. B 26, 101 (2002). 
[8] V. Apinyan and R. Melin, Eur. Phys. J. B 25, 373 (2002). 
[9] M.J.M. de Jong and C. W. Beenakker, Phys. Rev. Lett. 74, 1657 (1995). 
[10] R.J. Soulen et ah, Science 282, 85 (1998). 
[11] S.K. Upadhyay et al., Phys. Rev. Lett. 81, 3247 (1998). 
[12] P. Tedrow and R. Meservey, Phys. Rev. Lett. 26, 192 (1971); 

P. Tedrow and R. Meservey, Phys. Rev. B 7, 318 (1973); 

R. Meservey and P.M. Tedrow, Phys. Rep. 238, 173 (1994) and references therein. 
[13] P. Fulde and A. Ferrel, Phys. Rev. 135, A550 (1964). 
[14] A. Larkin and Y. Ovchinnikov, Sov. Phys. JETP 20, 762 (1965). 
[15] M.A. Clogston, Phys. Rev. Lett. 9, 266 (1962). 

[16] E.A. Dernier, G.B. Arnold and M.R. Beasley, Phys. Rev. B 55, 15174 (1997). 

[17] A.I. Buzdin, L.N. Bulaevskii, and S.V. Panyukov, JETP Lett. 35, 178 (1982) [Zh. Eksp. Teor. Fiz. 35, 147 (1982)]. 

A. Buzdin, B. Bujicic, and M. Yu. Kupriyanov, Sov. Phys. JETP 74, 124 (1992) [Zh. Eksp. Teor. Fiz. 101, 231 (1992)]. 
[18] V.V. Ryazanov, V.A. Oboznov, A. Yu. Rusanov, A.V. Veretennikov, A.A. Golubov, J. Aarts, Phys. Rev. Lett. 86, 2427 

(2001). 

[19] T. Kontos, M. Aprili, J. Lesueur, and X. Grison, Phys. Rev. Lett. 86, 304 (2001). 
[20] A.I. Buzdin and M. Yu. Kupriyanov, JETP Lett. 52, 487 (1990); 

A.I. Buzdin, M. Yu. Kupriyanov and B. Vujicic, Physica C 185 - 189, 2025 (1991). 
[21] J.S. Jiang, D. Davidovic, D.H. Reich, and C.L. Chien, Phys. Rev. Lett. 74, 314 (1995). 
[22] C.L. Chien, J.S. Jiang, J.Q. Xiao, D. Davidovic, and D.H. Reich, J. Appl. Phys. 81, 5358 (1997). 

[23] L.V. Mercaldo, C. Attanasio, C. Coccorese, L. Maritato, SL. Prischepa, and M. Salvato, Phys. Rev. B 53, 14 040 (1996). 
[24] J.S. Jiang, D. Davidovic, D.H. Reich, and C.L. Chien, Phys. Rev. B 54, 6119 (1996). 

[25] Th. Muhge, N.N. Garif'yanov, Yu. V. Goryunov, GG Khaliullin, L.R. Tagirov, K. Westerholt, LA. Garifullin, and H. 
Zabel, Phys. Rev. Lett. 77, 1857 (1996). 

Th. Muhge, K. Westerholt, H. Zabel, N.N. Garif'yanov, Yu. V. Goryunov, LA. Garifullin, and GG Khaliullin, Phys. Rev. 
B 55, 8945 (1997). 

[26] M.D. Lawrence and N. Giordano, J. Phys. Condens. Matter 39, L563 (1996). 

[27] V.A. Vas'ko, V.A. Larkin, PA. Kraus, K.R. Nikolaev, D.E. Grupp, CA. Nordman, and A.M. Goldman, Phys. Rev. Lett. 
78, 1134 (1997). 

[28] M. Giroud, H. Courtois, K. Hasselbach, D. Mailly, and B. Pannetier, Phys. Rev. B 58, R11872 (1998). 
[29] V.T. Petrashov, LA. Sosnon, I. Cox, A. Parsons, and C. Troadec, Phys. Rev. Lett. 83, 3281 (1999). 
[30] AT. Filip, B.H. Hoving, F.J. Jedema, B.J. van Wees, B. Dutta, and S. Borghs, Phys. Rev. B 62, 9996 (2000). 
[31] M. Giroud, K. Hasselbach, H. Courtois, D. Mailly and B. Pannetier, cond-mat/0204140. 
[32] I. Baladie, A. Buzdin, N. Ryzhanova, and A. Vedyayev, Phys. Rev. B 63, 054518 (2001). 
I. Baladie and A. Buzdin, cond-mat/0209466. 



20 



[33] P.G. de Gennes, Phys. Letters 23, 10 (1966). 

[34] G. Deutscher and F. Meunier, Phys. Rev. Lett. 22, 395 (1969). 

[35] J.J. Hauser, Phys. Rev. Lett. 23, 374 (1969). 

[36] E. Vecino, A. Martin- Rodero and A. Levy Yeyati, Phys. Rev. B 64, 184502 (2001). 

[37] P.G. de Gennes, Superconductivity of Metals and Alloys (Benjamin, New York, 1966). 

[38] N. Stefanakis, Phys. Rev. B 66, 024514 (2002). 

[39] P.G. de Gennes and D. Saint- James, Phys. Lett. 4, 151 (1963). 

[40] J.C. Cuevas, A. Martin-Rodero, A. Levy Yeyati, Phys. Rev. B 54 736 (1996). 

[41] K. Halterman and O.T. Vails, Phys. Rev. B 65, 014509 (2002). 

[42] M. Tinkham, Introduction to superconductivity, Second Edition, McGraw-Hill (1996). 

[43] G. Deutscher and F. Meunier, Phys. Rev. Lett. 22, 395 (1969). 

[44] J.J. Hauser, Phys. Rev. Lett. 23, 374 (1969). 

[45] Y.A. Genenko and Y.A. Ivanchenko, Theor. Math. Phys. 69, 1056 (1986). 

[46] N. Moussy, H. Courtois and B. Pannetier, Europhys. Lett. 55, 861 (2001). 



21 



